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ABSTRACT 


A review of the results published for systems containing 
hydrogen sulphide, carbon dioxide, methane, nitrogen, and 
natural gas is presented. The progress of the concept of 
vapour-liquid equilibrium is studied with a brief discussion 
on the methods of correlation and theoretical prediction of 
the vapour-liquid equilibrium data. 

The equilibrium cell was tested for its performance on 
multicomponent systems. Some modifications were incorporated 
and a change in the design of the piston is suggested. 

The two phase envelope of the quaternary natural gas 
system containing hydrogen sulphide, carbon dioxide, methane, 
and nitrogen was defined by measuring percent liquid present 
for different pressures at nine temperatures, 35°, 60°, 90°, 
ZO a laO? ESOS Gor, 170°, and 80°F over a range of 
pressures from 246 to 1,945 psia. The phenomenon of critical 
opalescence was observed. A comparison is made with the 
observations reported by other authors. A photograph showing 
critical opalescence is presented. 

Equilibrium ratios for hydrogen sulphide, carbon dioxide 
and methane were calculated at 120°F. These values were 
compared with those presented by three other authors on allied 
systems. The K-values for methane and carbon dioxide compare 
reasonably well, but larger discrepencies were observed for the 


hydrogen sulphide K-values. 
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A. INTRODUCTION 


A knowledge of the composition and amounts of the 
coexisting phases under heterogeneous equilibrium and the 
ability to determine the thermodynamic properties of fluids 
is basic for most convenient and economic use of the available 
resources. 


Souders, Selheimer, and Brown (166) | 


developed the 
concept of ‘equilibrium constants', and it has been found to 
be a suitable measure of the vapourization characteristics 
of petroleum mixtures (65) , 

They arbitrarily defined the equilibrium constant K, 
as the ratio of the mole fraction of a constituent in the 
vapour phase, y, to the mole fraction of that constituent in 
the liquid phase, x, at a defined temperature and pressure: 


of the equilibrium. Therefore, at equilibrium 


Kp cay yee, (1) 


Equilibrium constants are functions of temperature, 
pressure, and the overall composition of the system. The 
effect of compositions is very small at low pressures, but, 
above about 1,000 psia, the equilibrium constants depend 
considerably on the overall composition of the system (167) , 


The term constant is thus a misnomer, but it is used widely 
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(107) 


in the petroleum industry. Muskat has proposed the term 
"Equilibrium Ratio'. 

Non-hydrocarbons such as hydrogen sulphide, carbon 
dioxide and nitrogen are important components of the fluids 
produced from underground petroleum reservoirs. Most binary 
systems containing paraffin hydrocarbons and some ternary and 
more complex systems have been investigated. Methods for 
correlating and predicting such data are available. However, 
the information available on systems of hydrocarbons and non- 
hydrocarbons such as hydrogen sulphide and carbon dioxide is 
very scanty. These systems, in general, exhibit a highly 
non-ideal behaviour. 

The equilibrium cell used in this work was constructed 
at the University of Alberta in 1961. sandercock (158) has 
presented a detailed description of the cell. Hughes (9) 
used the cell to study the phase behaviour of the binary 
system hydrogen sulphide-butane. 

The purpose of this investigation was: 

(1) to assess the performance of the existing equipment with 
a multicomponent system and to modify it, if necessary. 
(2) to design a convenient unit for making multicomponent 


gas samples of any desired composition. 
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(3) 


(4) 


to employ the equipment in the study of the phase behaviour 
of a natural gas system at a series of suitable pressure 
and temperatures. 

to compare the equilibrium ratios thus obtained with 


values reported in the literature on related systems. 
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B. THEORY 


(1) Literature Review 


The vapour pressures of hydrogen sulphide were reported 


by west (285) , (121) 


Reamer, Sage, and Lacey presented data 

on its phase and volumetric behaviour. Data on several binary 
systems containing hydrogen sulphide have been reported. The 
phase and volumetric behaviour of the system methane - hydrogen 
sulphide was reported by Reamer, Sage, Lacey (122) , They 
studied the system at temperatures ranging from 40 to 340°F 


(78) studied 


and pressures up to 10,000 psia. Kohn and Kurata 
this system for the temperature range of -300 to 300°F and 
pressures up to 2,000 psia and thus extended the data to lower 
temperatures. They report the existence of type K singular 
points for this system where two phases are in critical identity 
in the presence of a third phase. This system did not exhibit 
any azeotropic behaviour for the whole range of temperatures 
investigated. 


Kay and Brice (7) 


reported the system ethane-hydrogen 
sulphide for the temperature range of -7 to 212°F and pressures 


up to 1,200 psia. A minimum boiling azeotrope was observed 


at concentrations high in ethane and persisted up to the 
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critical point. The phase equilibrium data for the system 


hydrogen sulphide - propane were investigated by Kay and 


Q 
i) 


“ 6 = eal E 
Rambosek (9 ) from 32°F to the critical locus. Brewer, Rodewald 
? n 1 ¥ 0 ~ 0 ’ 7 
and Kurata (32) extended the data on this system to the solid - 
liquid - vapour region which starts about -125°F. They 
studied this system for the temperature range of -130 to 
140°F and verified the data of previous authors in the region 
ae ee , (169) 
common to both the investigations. Steckel has studied 
this system in the temperature range of -30 to 15°C while 
. . hn J s (54) see: s e ° s A. 
Gilliland and Scheeline restricted their investigation to 
the temperature range of 124 to 201°F and pressures between 
490 and 600 psia. In this system, an azeotrope was formed 
at high hydrogen sulphide concentrations and persisted up to 
the critical region. 


(158) 


Sandercock studied the phase behaviour of the 


system hodrogen sulphide - butane at 100°F and Hughes (6°) 
has reported data in the temperature range of 125 to 275°F and 
pressures up to 1,780 psia. No actual azeotrope was observed 
but his data indicate a tendency toward an azeotropic behaviour. 


Reamer, Sage, and Lacey ‘123) 


studied the phase and 
volumetric behaviour of the system pentane - hydrogen sulphide 


in the temperature range of 40 - 460°F and for pressures up to 


10,000 psia. No azeotrope was formed. The phase and volumetric 
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behaviour of the system hydrogen sulphide - decane was reported 


2A. 
by Reamer, Selleck, Sage, and nacey (+24) 


at temperatures 
ranging from 40 - 460°F and for pressures up to 10,000 psia. 
They found the system to be of the common binary type with 
no unusual behaviour. 

The hydrogen sulphide - carbon dioxide binary has been 


(23,24) 


studied by Bierlein and Kay from 0°C to the critical 


point of hydrogen sulphide. The system exhibited a tendency 
to form a minimum boiling azeotrope, but no actual azeotrope 
‘ 4 (169) ee 2 

is reported. Steckel has reported isothermal dew and 
bubble point pressures for this system at 0°, -26.8°, and 
-52°C. His results at 0°C, though qualitatively similar 
differ somewhat quantitatively from those of Bierlein and 


(87) 


Kay. Kurata and Sobocinski investigated this system from 


the critical region to the solid-liquid-vapour region. 


140) 


Robinson et al (439, investigated the methane - 


carbon dioxide - hydrogen sulphide ternary. These investigators 
reported equilibrium constants for each component as a function 
of pressure and ratio of other components at 100°F, 50°F, 

and 160°F and at pressures up to 1,600 psia. 


(42) 


Very recently Cramer has compiled the phase 


diagrams for carbon dioxide up to 12,000 atmospheres. muenan (23-86) 


was perhaps the first to investigate the behaviour of paraffin 


feaoger sot ana 
Lciehiainas sqymiot 


fy 


ARAL a eaO a, ot ob ee ore 


‘tine “esc ali °0 ie mageye aut sot soma 


hydrocarbons and carbon dioxide. His measurements indicated 
a minimum boiling azeotropic behaviour for the system ethane - 
carbon dioxide. sweigert (172) has compiled the thermodynamic 
properties of carbon dioxide in the temperature range of -75 
to 1,800°F and for pressures up to 3,000 psia. 

Reamer, Olds, Sage and Lacey (125) investigated the 
system methane —- carbon dioxide in the gaseous region at 
temperatures ranging from 100° to 460°F and throughout the 
pressure interval from 100 to 10,000 psia. Donelly and 
Katz (47-48) have reported the solid-liquid-vapour equilibrium 
data for this system in the temperature range of -100 to 29°F 
and for pressures up to 1,150 psia. 

The volumetric behaviour for the system ethane - 
carbon dioxide has been reported by Reamer, Olds, Sage, and 


Lacey (426) 


at temperatures ranging from 100 - 460°F and for 
pressures up to 10,000 psia. 

Poettman and Katz (118) reviewed the literature relating 
to binary mixtures of carbon dioxide and hydrocarbons. In 
addition, they studied the phase behaviour of mixtures of 
carbon dioxide with propane, butane, and pentane. They 
investigated the system propane — carbon dioxide at temperatures 


ranging from 32 to 210°F and for pressures up to 1,000 psia, 


the system butane — carbon dioxide in the temperature range 
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of 45 to 300°F and for pressure up to 1,180 psia, and the 
system pentane - carbon dioxide at temperatures between 45 
and 365°F and for pressures up to 1,450 psia. 


Reamer, Sage, and Lacey (127) 


presented data on the 
volumetric behaviour of the system carbon dioxide - propane 
at temperatures ranging from 40 to 460°F and pressures up to 
10,000 psia. They studied the vapour-liquid equilibrium of 
this sytem in the temperature range of 40 to 160°F and for 
pressures up to 1,002 psia. Akers, Kelley, and Lipscomb ‘1? 
worked at lower temperatures between ~-40 and 32°F and for 
pressures up to 507 psia. The critical locus of this system 
shows a sharp downward curvature at the carbon dioxide end of 
the line, thereby indicating that the forces tending to form 
a minimum boiling azeotrope as in the ethane ~- carbon dioxide 
system have not yet completely disappeared. 

Olds, Reamer, Sage and ieieeuwe cae investigated 
volumetric behaviour of the system butane - carbon dioxide 
at temperatures between 100 and 460°F and for pressures up to 
10,000 psia. They also reported data on the vapour-liquid 
equilibrium of this system in the temperature interval 100 - 
280°F and for pressures up to 1,200 psia. This system also 
exhibited a reversal in the curvature of the bubble point 


pressure-composition relationships. 
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Mathews anda Hura(99) in 1946 utilized the data of 
Kvalnes and Gaday (88) and Olds, Reamer, Sage, and Lacey (110) 
to establish the P-V-T behaviour of methane in the superheated 
region. They also used the data of Keyes, Taylor, and smith (7°) 
on the saturated liquid densities of methane in conjunction 
with the Clapeyron equation to develop the densities for the 
saturated envelope. This enabled them to establish densities 
up to Pr = 20. The data of Schamp, Mason, Richardson, and 
Altman (161) | Michels and Nederbragt (104,105) | Keyes and Burke (74) | 
Keyes, Smith, and Joubart (75) , and those of Amagat (6. 7) reinforce 
the information resulting from the data of Olds, Reamer, Sage, 
and Lacey and Kvalnes and Gaddy in the intermediate to high 
pressure region. The data of Kazarnovskii and Levchenko (69) 
extended the high pressure data of Olds et al to temperatures 


(100) calculated the 


up to Tr = 3.08. Matschke and Thodos 
densities of compressed liquid methane up to pressures of Pr = 20 
and temperatures as low as Tr = 0.5. 

The system methane - ethane has been investigated by 


Ruhemann (147) 


in the temperature range from -104 to 0°C 
and pressures up to 100 atmospheres. At pressures below 15 


atmospheres, the system behaved nearly ideally obeying the 


perfect solution laws. 
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Sage, Lacey, and Schaafisma (149) 


investigated the system 
methane - propane at temperatures between 20 to 90°C and pressures 
up to 200 atmospheres. The system was found to have a large 
region of retrograde condensation of the first kind. Akers, 
Burns, and Fairchild 2) extended the data on the system to 
lower temperatures ranging between -176 to 32°F and for pressures 
up to 1,450 psia. 

The system butane - methane has been studied extensively 


(108,128,150,151,137) ame combined data 


by various authors 
of these investigators embraced the range of temperatures from 
-4° to 460°F and pressures from 147 to 10,000 psia. 

Sage, Reamer, Olds, and Lacey (153) have investigated 
the phase and volumetric behaviour of the system methane - 
pentane at temperatures ranging from 100 to 400°F and for 
pressures up to 5,000 psia. 

The methane - decane system has been investigated 
over a wide range of temperatures from 70 - 460°F and for 
pressures up to 10,000 psia (152,90,153) | 


Volumetric measurements on nitrogen have been reported 


for temperatures from 0 to SO00°F at pressures from atmospheric 


POs Moodmesia|( 13.1190), 2 14/59/8102, 103,62), (154) 


Sage and Lacey 
and Bloomer and Rao (30) have tabulated the thermodynamic prop- 


erties of nitrogen. 


AGIs 49tq 20% bas 1808 a £ ~ineiounted 


ylevianed#s batbwte aged s 


Akers, Attwell, and Bobineon have investigated the 
system nitrogen - butane. The system nitrogen - heptane was 
reported by Akers, Kehn, and Rlgeren . Some data are 
available for the system nitrogen - carbon dioxide at 0°C 
and for pressures up to 7,000 ae, O28 

Various investigators have reported data on systems 
containing a natural gas and hydrocarbon liquids such as 
crude oil, absorption oil, and gas distillate. In some 
cases, they have reported the equilibrium ratios for components 
such as hydrogen sulphide, carbon dioxide, methane, and 
nitrogen. 

The systems containing hydrogen sulphide, carbon 


dioxide, methane, nitrogen, and natural gas on which data 


are available are summarized in tables 1 to 5 respectively. 


Table No. 1 


REFERENCES ON HYDROGEN SULPHIDE SYSTEMS 


CoHe-Hos 


C3Hg-HoS 


Cai 9-H2s 


n-C5H12-H2S 
N-C 1 9Ho9-HoS 


H9S-CO5 


Ho S- H),-CO> 


Authors 


Reamer, Sage, and Lacey 
Kohn and Kurata 


Kay and Brice 

Kay and Rambosek 

Brewer, Rodewald, and Kurata 
Gilliland and Scheeline 


Steckel 


Hughes 
Sandercock 


Reamer, Sage, and Lacey 


122 
78 


67 


Reamer, Selleck, Sage, and Lacey 124 


Steckel 


Bierlein and Kay 
Kurata and Sobocinski 
Steckel 


Robinson, Lorenzo, and 
Macrygeorgos 

Robinson and Bailey 

Kennedy and Wieland 


Selleck, Carmichael, and Sage 


H,S-H,0-ethanolamine Atwood, Arnold, and Kindrick 


HoS-Reservoir fluid Vagtborg 


163 


181 
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Table No. 2 


REFERENCES ON CARBON DIOXIDE SYSTEMS 


system 


€1-Co2 


Co-CO5 


C3-CO> 


Cy-CO> 


C5-CO5 


H5S-CO> 


CHy-CO9-H2S 


C2H2-CO2 

C2Ha-Co2 

(C5Hs5) 20-COo 

CO5-SO5 

CO5-Natural Gas- 
Condensate 


CO5-Natural Gas- 
Crude Oil 


Authors 


Donnelly and Katz 
Reamer, Olds, Sage, and Lacey 
Tsikilis 


Reamer, Olds, Sage, and Lacey 
Kuenan 


Akers, Kelley, and Lipscomb 
Reamer, Sage, and Lacey 
Poettman and Katz 


Poettman and Katz 
Olds, Reamer, Sage, and Lacey 


Poettman and Katz 


Bierlein and Kay 
Kurata and Sobocinski 
Steckel 

Thiel and Schulte 


Citation 


47, 48 
125 
178 


126 
83) = 86 


Robinson, Lorenzo, and Macrygeorgos 140 


Robinson and Bailey 
Kennedy and Wieland 


Kuenan 

Tsikilis 

Thiel and Schulte 
Caubet 

Thiel and Schulte 


Poettmann and Katz 


Poettman 


126 
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CO5-N5 


C5H5C1-CO5 
€O5-Cl'5 
CO5- (CH3) 20 


CO9-air 


CO>-H> 
CO5-HC1 


CO2-CH3C1 


CO2-N20 
CO>-CgHeNO> 


CO2-02 


Table No. 2 (Continued) 


Tsikilis 
Torocheshnikov 
Andrews 

Thiel and Schulte 
Thiel and Schulte 
Thiel and Schulte 


Cailletet 
Van der Waals 


Vershaffelt 
Ansdell 


Kuenan 
Caubet 


Caubet 


Kohnstamm and Reeders 


Keesom 


78 
£75 

10 
174 
174 
174 
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182 


183 
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Table No. 3 


REFERENCES FOR METHANE SYSTEMS 


system 


CH4- CoH, 


CHy4-C3Hg 


CHg-n-C4Hi9 


CHy-i-C4Hio 


CHa-n-CsHy92 


CHy-i-CsH19 
CHy~hexanes 


CHy-nC7Hig 


CHy-nCi9H22 


CHa- CoHe-CoH, 
CHsa-CoHa 
CH4-C3H,-nC4Hi9 


CH4a-nC4aHji9-nCj9H22 


Authors 


Ruhemann 


Bloomer, Gami, and Parent 
and Schaafasma 
and Lacey 


and Fairchild 


Sage, Lacey, 
Reamer, Sage, 
Akers, Burns, 


Neberbraght 

Reamer, Korpi, Sage, 
Sage, Budenholzer, and Lacey 
Sage, Hicks, and Lacey 
Rigas, Mason, and Thodos 
Olds, Sage, and Lacey 
Sage, Reamer, Olds, 
Taylor, Wald, Sage, 
Boomer, and Johnson 


and Lacey 
and Lacey 


Amick, Johnson, and Dodge 
Boomer and Johnson 


Boomer and Johnson 
Reamer, Sage, and Lacey 


Sage, Lavender, and Lacey 
Reamer, Olds, Sage, and Lacey 


Guter, Newitt, and Ruhemann 
Guter, Newitt, and Ruhemann 
Rigas, Mason, and Thodos 


Reamer, Fiskin, and Sage 
Reamer, Sage, and Lacey 


and Lacey 


Citation 


e29 


B52 
153 


56 


56 
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130 
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14. CHy~nC,Hy ~H,0 McKetta and Katz Low 
15. CHy-H20 Olds, Sage, and Lacey ear2 
16. CHg-No9-Cyclohexane Johnson 64 
17. CH4a-N2-Heptane Johnson 64 
18. CH4g-N9-Benzene Johnson 64 
19. CHyg-No-Toluene Johnson 64 
20. CHyg-Crystal Oil Sage, Backus, and Lacey BSS 
dk CH,-Kensol-16 Rzasa 148 
22. CHy-COs Donelly and Katz 47, 48 
Reamer, Olds, Sage, and Lacey 25 
Tsikilis 178 
23. CH,-H 5S Reamer, Sage, and Lacey u22, 
Kohn and Kurata 78 
24, CHyg-CO5-H>S Robinson, Lorenzo, and Macrygeorgos 140 
Robinson and Bailey 139 
Kennedy and Wieland US 
2% CH,-NH3 Krichevsky and Tsikilis 80 
26%. CH,-CO Toyama 176 
An. CHi~—Crade Oi1= 
Natural Gas Buckley 34 
Zoi CH,-Natural Gas Davis, Bertuzzi, Gore, and Kurata 45 


29. CHyg-Crude Oil Sage, Webster, and Lacey 156 
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Table No. 4 


REFERENCES ON NITROGEN SYSTEMS 


system 
No-nCgHyi9 
No-nC 716 


Nog-CO5 


No—CO9-O09 
No-CO 
N2-C2Ha4 
N5-H 0 
No-SO5 

No- Benzene 


No-NH3 


No-Natural Gas 
No-He 


No-A 


Authors 
Akers, Attwell, and Robinson 
Akers, Kehn, and Kilgore 
Tsikilis 
Torocheshnikov 
Andrews 
Zenner and Dana 
Schiller and Canjar 
Tsikilis 
Bassett and Dodge 
Tsikilis 
Krishchevsky 
Wiebe and Gaddy 
Krishchevsky and Tsikilis 


Lindroos and Dodge 
Krishchevsky and Bolshakov 


Davis, Bertuzzi, Gore, and Kurata 


Buzyna, Macriss, and Ellington 


Long and DiPaolo 
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Table No. 5 


REFERENCES ON NATURAL GAS SYSTEMS 


system 


Natural gas 


Natural gas- 
distillate 


HoS-COo-Natural 
Gas 


neC7h,¢ ~Natural 
Gas 


Natural gas-Np or 
CH, 


COo-Natural gas- 
Crude Oil 


CHy-Natural Gas- 
Crude Oil 


Natural gas-Crude 


Oil 


Natural gas- 
Absorption Oil 


Natural Gas-Natural 


Gasoline 


Multicomponent 
hydrocarbon 
systems 


Hydrocarbons in 
Absorption Oil 


Authors 
Gore, Davis, and Kurata 


Brown and Stutzman 


Poettman and Katz 
Roland, Smith, and Kaveler 


Robinson, Macrygeorgos and Govier 


Boomer, Johnson, and Piercey 
Davis, Bertuzzi, Gore and Kurata 
Poettman 

Roland 

Buckley 

Katz and Hachmuth 

Standing and Katz 

Jacoby and Rzasa 

Kirkbride and Bertetti 
Jacoby and Rzasa 

Katz and Kurata 


Sage and Lacey 


Webber 


Citation 


315) 
33 


116 
145 


141 


cl 


45 
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144 


168 


184 
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(2) Phase Behavior Principles 


Equilibrium in a close system implies a situation in which 
there is no change with respect to time. If we consider a 2-phase 
system at equilibrium containing n components and concern 
ourselves with the mole fraction y; of the gee component in 
the vapour phase and x; of the same component in the liquid 
phase, the equilibrium ratio for this component is arbitrarily 
defined as 

Ki = y4/%3 “ 

The definition of "ideal equilibrium ratios" comes from 
the concept of a "perfect solution". If the molecules of the 
components of the solution are so similar that the forces 
between unlike molecules are the same as those between like 
molecules and if the components mix without the complicating 
effects of molecular association or chemical combination the 


solution would obey Raoult's eagle), 


where, 


th 


Puls che pantsal pressure of the i component over 


the solution, and 


Py i is the vapaur pressure of the jth component in its 


pure state at the same temperature. 
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Dalton's Law (44) states that the total pressure of a 
gas is equal to the sum of the partial pressures of the com- 


ponents present, thus 
Py Ya (3) 


where, 
P is the total pressure of the system. 
Equations (2) and (3) may be combined to give, 


Pyi 


K- = V4/%s = é 


ab 


This equation suggests a useful method of representing 
the equilibrium ratios graphically. Taking logarithms of both 


sides of equation (4), one obtains 


lege) > 10S) Py. =) deg, P (5) 


Hence a log-log plot of K; versus P will give a straight 
line having a slope of -l, passing through K; = 1 where P = Pv 
at the temperature in question. 

It has been shown that Raoult's and Dalton's Laws are 
not exact even in relatively simple mixtures of hydrocarbons 
(37,142,192) - 

In general, K—factors predicted by this method are 
restricted to low pressures, some considerable distance from 


the critical pressure and to ideal solutions. Normal paraffin 
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hydrocarbons form ideal solutions over certain ranges of 


pressure and temperature (143) | 


but the presence of components 
such as water, hydrogen, carbon dioxide, hydrogen sulphide, 
or olefins frequently cause significant deviations. 

G.N. Lewis (95) proposed the use of fugacity in place 
of pressure whenever the fluids behave non-ideally. A closer 
approximation to actual k-factors is obtained if partial 
pressure and pressure are replaced by fugacities in making 
equilibrium calculations. An ideal solution is defined as 
one in which the fugacity of each component is proportional 
to the mole fraction of that component at every temperature 


and pressure (94) , 


The equilibrium ratios defined in terms of 
fugacity are sometimes called ‘ideal solution equilibrium 


Gatios'. 
Defining the fugacity at constant temperature T, as 
ain eo =" "v dp (6) 
where, 
f is the fugacity, v is the molal volume of vapour, 


and p the partial pressure of the component. The Dalton's 


Law type equation becomes, 


Ee enV a (7) 


and the Raoult law type equation becomes 
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where, a bar denotes the fugacity of component in the mixture, 
v and L stand for vapour and liquid respectively. 


The criteria for equilibriun, which are 


V L V L ame NY / Sous 
T SOoPu sp = p , and 54, = 4 (9) 
give 
L 
ndeX.; = 
Ki deal = 24 = = (10) 
eal, al 


In this case, fugacities have replaced the pressure in 
equation (4). K. geal is a function of system temperature 
and pressure only and is independent of composition. 

Many generalized K-factor correlations based on the 


4 fe) 
g (20,49, 119) | These 


fugacity concept have been presente 
idealized approaches are not truly descriptive and large 
deviations are encountered in non-ideal systems. Experimental 
data have shown that the K-factors for any component are 
functions of the temperature, the pressure, and the composition 
of the system. Various attempts have been made to devise a 
suitable parameter to account for the effect of composition. 
The K-factors are observed to converge to unity at some 
pressure for each temperature chosen. This leads to the 
concept of convergence pressure, which at the present time is 


generally accepted as the best means available to relate 


K-factors to composition. 
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The convergence pressure is defined as the pressure 
at which K's for all components appear to converge to unity 
at the system temperature. The K's are actually continuous 
up to K=1l only when the equilibrium temperature happens to 
be the critical temperature of the system. At temperatures 
other than the critical, the discontinuous curves on a log K 
vs log P plot when extended appear to converge to Gindey ee 

From Gibbs phase rule, for two phase equilibria, 

K = @ (P,T, Composition of (n-2) components) (11) 
where, n is the number of components of the mixture. 


If we express, 


Wisriei- APs Bee) (12) 
then 
Poy = 6" (T, Composition of (n-2) Components) (13) 
where, P., is the convergence pressure of the system. For 


binary systems the convergence pressure is the critical 
pressure of the mixture having a critical temperature equal 
to the system temperature. 

The application of convergence pressure to multi- 
component systems requires treating the multicomponent mixture 
as a hypothetical binary while estimating the K-factors. 


Many methods have been presented to calculate the con- 


(57,95,165) (58,193) 


vergence pressure and charts are available 
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to obtain the equilibrium ratios when the convergence pressure 
is known. 

For mixtures of paraffin hydrocarbons the convergence 
pressure is found to be a satisfactory parameter to account 
for the effect of composition. However, the addition of non- 
hydrocarbons may destroy the usefulness of the convergence 
pressure type correlations. 

Activity coefficients have been used in the case of 
reasonably ideal systems to extend the vapour-liquid equilibrium 


(52,164) Smith and Waceon one employed 


data to high pressures 
high pressure activity coefficients to correlate the vapour- 
liguid equilibrium data on paraffin hydrocarbons. In the 184 
systems they studied an average deviation of 15% was observed. 
The mean deviation was 216% without the use of activity 
coefficients. 

Normal boiling points have been used as a correlating 
parameter for the vapour-liguid equilibrium data on paraffin 
hydrocarbon eestenee but their application is limited 
only to such systems. 

The techniques available to predict equilibrium ratios 
are 


(i) Raoult's and Dalton's Laws. Their application is 


limited to systems approximating ideal solutions. 
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(ii) Cluster theory of vapour-liquid equilibria. Lu, 
(22) 


James, and Ting postulated that deviations from the ideal 
behaviour are due to the formation of clusters in liquid and 
vapour phases. 

By writing equilibrium equations between clusters and 
molecules of the two components, these authors calculated the 
vapour-liquid equilibria for binary systems. Assuming the 
clusters to be formed among three components or more, the 
above method was extended to calculate the vapour-liquid 
equilibria for complex systems. Applicability of this theory 
has not been tested sufficiently so far, particularly so for 
hydrocarbon - non hydrocarbon systems. 

(iii) Activity coefficients. One can define the deviation 
factors from Raoult's law as follows: 


Pp . 
is peaeeeins (14) 


Py. “ia 
The term introduced as a deviation factor from Raoult's 
law can be shown to be equivalent to the thermodynamic 
property, activity coefficient. 
Lewis and Randall (94) defined activity, as the relative 
fugacity, or the ratio between fugacity of the substance in 
solution and its fugacity at a standard state. Thus, for pure 


liquid as the standard state; 
fi 
£4 


(15) 
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Activity coefficients are simply the activities 


divided by their respective mole fractions of ¥, 
i 


Salt , TAS) s 
x j 


i= — (16) 


eee x4 


Where the vapours are perfect gases, £; = p; and 


a Pyss so that 


It is thus apparent that the deviation factors defined 
by equation (14) are equivalent to values of activity coefficients 
under conditions that the vapours approximate perfect gases. 


A number of methods have been proposed for predicting 


rs ee 2 (aS 3Si89), os), 135, 136,159 7194) ‘ 
the activity coefficients with 


the objective to permit the evaluation of YY, as a function 
of liguid phase composition from a minimum of experimental 
data. These methods are derived from the solutions of the 
Gibbs'-Duhem equation when expressed in terms of activity 


coefficients. The most well-known of these solutions are the 


(89) (98) 


equations derived by van Laar Margules , and Scatchard 


and Hamer (159) , 


These may be used to calculate vapour-liquid equilibria 
for non-ideal systems. Only a knowledge of van Laar constants 


is needed. These equations have been employed in the case of 


i i i i i ’ ' gD 
non-associating solutions with considerable eHeeaBeie 41,51,53 
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For solutions of associating components, the van Laar 
equations are generally unsatisfactory. Similarly, with two 
coefficients the Margules relations are inadequate for most 
associating systems and particularly binary systems in 
which activity coefficients for the two components are non- 
symmetrical. In such cases, the van Laar type equation is 
usually superior. Margules type equations give better results 
when more terms are included. Scatchard and Hamer equations 
have advantages over others in some special cases. They fail, 
sometimes, due to the lack of accurate data on boiling points 
as a function of composition at constant pressure. These 
equations are valid only for constant temperature data. They 
have not been used for systems containing non-hydrocarbon 


components such as HS or CO.- 


(iv) Theory of liquid mixtures. Sweeny and Rose (171) have 
developed a very complex method of predicting vapour-liquid 
equilibria using a theory of liquid mixtures. Its applicability 
to hydrocarbon systems has not been studied. 


(v) Activity coefficients from molecular structure. Deal, 


46,115,192 
et ai ) 


suggested that in any solvent environment, 
the partial molal excess free energy of solution of a solute 


molecule can be taken as the sum of individual contributions 


from each of its structural groupings and that these contri- 
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butions depend in some way upon the numbers and kinds of 
structural groupings which make up the environment. No exten- 
Sion of a group approach to excess free energies or activity 
coefficients has, however, appeared recently. 

(vi) The law of correspanding states. The theory of corres- 
ponding states was suggested by van der Waals. Its application 
to multicomponent systems was based on the Lewis-Randall ideal 
solution rule approximation. The K-factors predicted by this 
method were thus greatly in error near the critical point of 


mixtures. Leland et aaa) 


have recently derived a rigourous 
definition of pseudocritical temperature and pressure. This 
new definition is valid for all vapour-liquid compositions 
and for non-polar substances, even in the critical region of 
: (176) ' 
mixture. Toyama employed this method to the CH,-CO 


system. The overall average absolute deviations were 11.98% 


EO iCO) and s92:795 tors CHy 


(vii) Modified van Laar equations. Recently, Black (70727728) 


has developed a new equation of state. Using this equation 

of state, he has modified the van Laar type of equations. The 
effect of association and inter-association between the 
molecules has been accounted for. The vapour-liquid equilibrium 
of complex systems can be predicted from binary coefficients 
alone. The equation, however, has not been applied to hydro- 


carbon - non hydrocarbon systems. 
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(viii) Equation of State. An empirical equation of state 


proposed by Benedict, Webb, and Rubin 


(15,20) is frequently 


used to predict the vapour-liquid equilibria. The use of B-W-R 


equation is based upon the assumptions that 


ae 


it is applicable to each of the components, to their 
mixtures, and to both the vapour and the liquid 
phases. 

that the eight constants of the equation are known 
for the pure components and that the constants for 
the mixture can be calculated from these values, 
using the recommended mixture rules. This equation 
gives satisfactory results for most light hydrocarbon 


mixtures (43,162,170,106,176) | Results of Toyama (176) 


indicate that the agreement between the actual and 
predicted K-factors is poor for hydrocarbon - non 


hydrocarbon mixtures. 
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C. EXPERIMENTAL STUDIES 


1) Experimental Apparatus 


A very complete and detailed description of the vapour- 
liquid equilibrium cell was presented by Sandercock (158) | 
Hughes (9) has reported a few changes in the cell assembly. 
Basically, it is a windowed cell of variable volume. Some 
of the outstanding features of this apparatus are; 

(1) Position of the vapour-liguid interface can be 
changed without a change in sample volume. This is accomplished 
by connecting the top and bottom of the cell respectively to 
two ends of a double acting pump. This enables removal of 
mercury from one end of the cell at the same rate that it is 
added to the other. There is a movable piston at the top of 
the cell and the sample contained between the piston and the 
mercury at the bottom can be moved to any desired level in 
the cell. 

(2) The bottom of the cell is also connected to a 
positive displacement pump and the volume of the sample in the 
cell can be increased or decreased at will by the withdrawal 
or addition of mercury to the cell. 

(3) The window of the cell is lighted and it can be 


viewed through a telescope by means of a periscope arrangement. 
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A pin protruding from the bottom of the piston serves as the 
reference level for all volumetric measurements. The dis- 
placement pump and the double acting pump have been calibrated 
for volume previously (198) . The expected accuracy of volu- 
metric measurements is believed to be 10.005 ml. 

(4) There is a small port below the window of the 
cell. Mercury or liquid samples can be sprayed through this 
port into the gas phase to facilitate the attainment of 
equilibrium. Some difficulty was encountered in spraying when 
very small liquid samples were present. In this case, considerable 
time was required to allow the liquid to drain down the cell 
after spraying. 

(5) Samples are taken from a port at the window level 
through a calibrated sample pump. 

The schematic arrangement of the cell and its accessories 
is shown in figure 1 and figure 2 shows a section through 
the cell assembly. 

A change was incorporated in the positioning of the 
sample pump. The sample pump was placed horizontally in the 
work prior to this investigation. In that position, it was 
possible for a slug of the sample fluid to be trapped in the 


pump. In order to avoid any contamination of the samples 


due to such a possibility it was decided to place the sample 
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ASSEMBLY VIEW OF EQUILIBRIUM CELL 


FIGURE 2 
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pump vertically. A new pump was constructed which has a 
micrometer attachment for an accurate measurement of the volume 
of the samples withdrawn. The detailed drawing of the sample 
pump is presented in figure 3. 

Unlike the previous arrangement (60,158) | only one port 
is provided for the introduction of any fluid in the cell and 
also for the withdrawal of the samples. This eliminates some 
unnecessary lines. 

The temperature of the cell was measured by means of 
an iron-constantan thermocouple and a Leeds and Northrup 
model K-2 potentiometer. The thermocouple calibration is 
presented in figure 13 of the appendix. The expected 
accuracy of temperature measurements is *0.05°F. The tempera- 
ture of the reference junction was 32°F. 

The pressure in the cell was measured by means of two 
Heise bourdon tube gauges. The rating of these gauges were 
1,000 and 5,000 psi, and they were marked at intervals of 2 
and 10 psi respectively. The gauges were calibrated by means 
of a dead weight tester. The 1,000 psi gauge showed no error 
between the actual pressure and the scale reading for the 
entire range. In the range of pressures up to 3,000 psi, the 
maximum error in the scale reading of the 5,000 psi gauge was 


less than 12 psi. A constant correction was applied for the 
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mercury head from the gauge zero level to the cell level. 

The bath liquid used was ethylene glycol. The tempera- 
ture was controlled using an Aminco bimetallic type thermo 
regulator which was connected to a heater and refrigerator 


assembly by means of an automatic relay. 


(2) Experimental Procedures 


(a) Cell _ assembly:- Prior to installing the equipment 
for this investigation, the cell was dismantled and cleaned 
thoroughly. Some mercuric sulphide had settled on the window 
of the cell and great difficulty was experienced in removing 
it. In view of the difficulties that had been experienced 


(738) in the cell assembly, it was not advisable 


previously 
to take the windows out of the cell body for cleaning. After 
many unsuccessful attempts with several solvents it was found 
that an equimolal mixture of ammonium hydroxide and ammonium 
sulphide was an effective cleanser for removing the deposits. 
(b) Preparation of experimental mixture:- Recently a 


natural gas well in the Province of Alberta has been discovered 


which has the following composition: 


% by Volume 


Hydrogen Sulphide 68.00 
Carbon Dioxide 9...36 
Methane 21.64 
Nitrogen 1.00 
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It was decided to study the phase behaviour of a 
quaternary having the above composition. Mixing the four 
gases in a proportion close to the above provided great 
difficulty. Initially, an attempt was made to mix the four 
gases in a separate windowed cell, but great difficulty was 
experienced during its transfer to the main cell. It was 
therefore decided to design a convenient mobile unit to make 
the samples of a desired composition. 

A 300 ml. stainless steel bomb rated at 1,000 psi 
was used. The top of this bomb could be connected to the 
sample port of the main cell assembly. The bottom of the 
bomb was connected through a T to a Ruska displacement pump 
of 250 ml. capacity and to two driers in series. The driers 
were filled with anhydrous calcium chloride and Drierite. 
The calcium chloride drier was connected to the cylinder of 
the gas that was to be injected in the cell. A Heise pressure 
gauge rated up to 1,000 psi was used to measure the pressure 
in the bomb. The bomb was immersed in a constant temperature 
bath. A mercury thermometer was used to measure the tempera- 
tures. Mercury was used as the confining liquid. The schematic 
arrangement of the unit is shown in figure 4. 

All gases were taken in at 100°F. The volume of each 


gas required to make the desired composition and the corres- 
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ponding pressure were previously calculated. Initially, the 
cell, the bomb, and the adjoining lines were evacuated and 
filled with mercury. The valve connecting the cell assembly 
to the bomb was closed. The driers were purged with hydrogen 
sulphide gas. The desired volume of the gas was sucked in the 
bomb at the known pressure and temperature. The bomb was 
then isolated from the drier unit and the gas in the bomb was 
pushed into the cell by actuating the mercury displacement pump. 
In the same way, carbon dioxide, methane, and nitrogen were 
introduced in this order. Hydrogen sulphide was introduced 
first because it was the major component and because it has 
the lowest vapour pressure of all the components. Carbon dioxide, 
which has the next higher vapour pressure over the H5S was 
immediately followed to avoid its liquifaction in the lines 
at higher pressures. Nitrogen, since it was the trace com- 
ponent was introduced last. 

After the hydrogen sulphide gas had been admitted 
great caution was excercised not to add an excess of any other 
component. Minor deficiencies in the amounts of carbon 
dioxide, methane, and nitrogen were made up in this order. 
Intermediate analyses of the cell mixture when it existed in 
a Single phase were carried out to determine how far the 
composition was from the desired. The method of analysis is 


presented later. 
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An attempt was made to use the available compressibility 
factor charts to calculate the amounts required and the corres- 
ponding conditions of each gas to make the desired mixture. 
The results of these calculations were not found to be correct 
and thus trial and error was the only method available to mix 
the gases in the right proportion. 

Once the mixture of a composition close to the desired 
had been made, the valve connecting the cell to the sample 
pump was firmly closed. The data were taken in two parts. 
Firstly, only the two phase envelope of the system was deter- 
mined by measuring the amount of liquid present at the known 
pressure and temperature of the equilibrium. Regular checks 
were made to see if any appreciable change had occurred in 
the composition of the cell mixture. The composition was 
found to remain constant all through this part of the investi- 
gation. The second part consisted in analysing the two phases 
at equilibrium to calculate the equilibrium ratios for each 
component. Data obtained in this part of the investigation 
provided a check on the data of part one. 

The temperature was recorded and maintained constant 
all the time as equilibrium was being attained by spraying 
the liquid and mercury through the bottom valve. This 


could be taken as the point at which further spraying of 
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mercury caused no subsequent change in the pressure. The 
equilibrium pressure was recorded. The volumes of liquid and 
gas present were measured by noting the volumes indicated on 
the double acting pump when mercury, liquid, and the piston 
pin end were on the reference level in the window. This pro- 
cedure was repeated at other pressures and temperatures. 

Samples of gas and liquid were removed through the 
sample pump and carefully expanded to atmospheric conditions 
in a gas-burette for chromatographic analysis. The size of 
the liquid to gas samples withdrawn was in proportion to the 
amounts present in the cell. This helped in keeping the com- 
position of the cell mixture reasonably constant. A slight 
change did occur in the overall composition of the cell mixture 
but this would not affect the K-values appreciably. The 
sampling procedure has been reported previously (69) , 

As the sample was being removed from the cell, mercury 
was added to it at the same volumetric rate. Pressure was 
thus maintained constant and the equilibrium was not 
disturbed. 

The dew point was noted by the first appearance of 
liquid. This was further checked by reversing the procedure 
and noting the disappearance of the liquid. Extrapolation of 


the plot of percent liquid versus pressure gave the bubble 
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point pressure. Near the critical conditions neither the 
dew point nor the bubble point was apparent, and much longer 


time was required to achieve equilibrium. 


(3) Materials 


Mercury was used as the displacement liquid. Since 
the combination of mercury and hydrogen sulphide, in the presence 
of moisture caused almost instantaneous blackening of mercury, 
considerable care was taken to insure that pure mercury was 
used. The procedure for the purification of mercury has been 
reported previously (158) , 

All gases were dried before introducing them into «= 
cell. The gases used were all over 99% pure as indicated 
upon a chromatographic analysis and no further purification 
was necessary. Carbon dioxide and nitrogen were cbtained from 
Liquid Carbonic. Matheson supplied the C.P. grade hydrogen 


sulphide and pure methane was obtained from the Phillips 


Petroleum Company. 


(4) Analysis 


The vapour and liquid samples were analysed using a 
Burrel Ky Kromo-tog. Two columns were used. A column packed 


with ansul ether on 25% silocel separated methane, nitrogen, 
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and air as one component, carbon dioxide, and hydrogen 
sulphide. A molecular sieve type 13-X column separated 
oxygen, nitrogen and methane. Molecular sieves block carbon 
dioxide and hydrogen sulphide and are temporarily degenerated 
upon prolonged exposure to these gases. A constant check on 
the calibration revealed that it had not changed due to the 
slight amounts of the acid gases introduced with each sample. 
Satisfactory separation was obtained for all components. 
Pure gases were used for the calibrations. By using 7 samples 
of different volumes, sample volume could be plotted against 
the peak heights, satisfactory linear curves were obtained for 
all components and for both the columns. Four known synthetic 
mixtures of all the four gases were made in a gas burette and 
analysed and it was found that for both the columns there is 
a negligible interaction among the various components. The 
volumes read from the calibrations always added up to within 
1% of the actual sampler volume. The chromatographic cali- 
brations are presented in figures 14 and 15 of the appendix. 


The column conditions were: 


I Column Ansul ether on 25% silocel 
Column length 7 feet 8 inches 
Carrier gas Helium at 79.3 ml/min. 


Detector Current 200 mA 
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Column Temperature 


Sample Size 


Eel Column 
Column Size 
Carrier Gas 
Detector Current 
Column Temperature 


Sample Size 


Molecular sieve, type 13X 
8 feet 4 inches 
Helium at 67 ml/min 
200 mA 
LOE S 


1.85 ml 


The retention times were as follows: 


ah Column: 
Methane, air, and nitrogen O27 5 Min. 
Carbon dioxide 125) man. 
Hydrogen sulphide 3.10 min. 
OA Column: 
Oxygen i O25 man. 
Nitrogen 26575) Min. 


Methane 


4.500 min. 
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D. EXPERIMENTAL RESULTS 


The composition of the system used for this investi- 
gation was: 


% by Volume 


Hydrogen sulphide 68.04 
Carbon dioxide 9.30 
Methane 21.470 
Nitrogen 0.96 

100.00 


The two phase envelope of the system was determined 
at nine wcemperatunes 35°, 60iSe 90°, 120°, 1409, 150°, 160°, 
170°, and 180°F. Figures 5 and 6 are plots of the data on 
pressure versus percent liquid coordinates and with temperature 
as the parameter. The nine temperatures define the phase 
envelope adequately. Figure 7 is the two phase envelope for 
the system on pressure-temperature coordinates showing the 
lines of constant percent liquid. Figure 7 was obtained by 
cross plotting the data of Figures 5 and 6. The results 


are tabulated in tables 6 - 9 of the appendix. 


Critical Poant 
The phenomenon of critical opalescence was observed 
at 140°F and 1,761 psia. Figure 8 is a photograph of the critical 


phenomenon taken through a periscope-telescope arrangement. 
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FIG. 5 ISOTHERMS FOR THE H,S-CO, - CH, -No 
SYSTEM 
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FIGURE 8 THE CRITICAL OPALESCENCE 


The critical phenomena exhibited the following 
characteristics: 

(1) A reddish-orange opalescence was observed with 
transmitted light. 

(2) Reflected light produced a bluish-white colour. 

(3) The thickness of the interface between liquid and 
vapour increased as the critical point was approached by 
changing the pressure at the critical temperature. The inter- 
face was undefined at the critical point. 

(4) The opalescence was more intense in the liquid 
phase when the pressure was reduced from a higher value to 
the critical pressure. When the critical point was reached 
by slowly increasing the pressure, the vapour phase exhibited 
a more intense opalescence. 

(5) The opalescence was not observed in the main body 
of the cell mixture. 

(6) The fluid in the immediate vicinity of the piston 
pin showed a more intense opalescence. Sometimes a very fine 
spray of liquid rose in a jet toward the piston pin. 

(7) The critical point was reproducible within = O25 
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The analysis of the phases was made at only one tempera- 
ture 120°F. The equilibrium ratios were calculated for methane, 
carbon dioxide, and hydrogen sulphide. The values for nitro- 
gen were not considered adequately reliable due to the very 
low proportion of nitrogen present in both the phases. The 
results are plotted in figure 9 and are tabulated in tables 
8 and 9 of the appendix. 

The composition of the cell mixture had changed to the 
following due to the withdrawal of about 24 samples of vapour 


and liquid from the cell. 


%_by Volume 


Hydrogen sulphide 59.18 
Carbon dioxide 225 
Methane 29766 
Nitrogen AE 

100.00 


The experimental data were consistent and reproducible. 
The k-values obtained during this investigation are 
compared with the following: 
(lL) Data of Kohn and Kurata (78) at 140°F and) for a 
= 1,760 for the system methane hydrogen sulphide. (Fig. 10) 


Pov 


(2) Data of Robinson et a1 (439) for the system methane - 
carbon dioxide - hydrogen sulphide at 40°F and for Cy = 0.14, 
P = 1,760, psia. § (Fag. il) 
(3) wesma (29) values for Po = 1,760 and at 120°FR. 
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BR. DISCUSSION 


(1) Experimental Apparatus 

The equilibrium cell was used for a multi-component 
system for the first time. It was found to be satisfactory 
for the purpose. 

The vertical positioning of the sample pump offered 
some advantages. It is adequate to have just one port for 
injecting any component in the cell and also for the sampling 
of phases from the cell. 

A change in the piston design is suggested after this 
work. In the present design, the total volume of the cell 
mixture cannot be viewed through the window. The fluid trapped 
between the bottom of piston (in its lowest position) and the 
top of the window cannot be seen. It was thus not possible 
to make any volumetric measurements when very high percentage 
(above 90%) of liquid was present. It is therefore suggested 
that instead of the 1/8" diameter pin constructed at the bottom 
of the piston, a rod of 3/8" 0.d. (diameter of the passage 
connecting the two sections of the cell) be built in the 
piston. The length of this rod should be such that it can 


be seen near the top of the window when the piston is in 
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its lowest position. The rod can then be tapered into a pin 
that ends at the centre of the window. This will avoid the 
possibility of liquid drops of any appreciable volume sticking 


to the tip of the piston. 


(2) Discussion of Results 


(a) Two phase envelope:- A large retrograde condensa- 
tion region is apparent from the phase diagram for the system 
shown in Figure 7. The extrapolated constant percent liquid 
lines of Figure 7 also appear to converge at the observed 
critical point at 140°F and 1,761 psi. 

Since it was not possible to measure the volumes of 
phases correctly when more than 90% liquid was present, a 
certain amount of uncertainty exists at 60°, 90°, and 120°F 
for lines of percent liquid higher than 90%. This region is 
shown dotted on the phase envelope of Figure 7. 

(ob) Critical point:- The critical phenomena observed 
during the course of this investigation was similar to the 
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(177) and Young \! 7”) . 


one reported earlier by Travers and Usher 
Their finding that "the opalescence is confined to that phase 


which is decreasing in volume through movement of the dividing 
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surface, or at least is most intense in that phase" is 
confirmed. 

There may have been some temperature gradients present 
between the piston pin and the main body of the fluid. This 
might be responsible for the unsteady fountain-like stream of 
droplets striking the piston point as seen in figure 8. 

It was not possible to observe the critical opalescence 
throughout the mixture because the two barrels of the double 
acting pump are not of exactly the same size; the little change 
of pressure that occurred when the position of mixture in the 
cell was changed disturbed the critical conditions appreciably. 

(c) 


Equilibrium ratios:- The equilibrium ratios were 


determined at one temperature only. The composition of the 
mixture had changed slightly due to the sampling of phases, 
but, it should not affect the K-values appreciably. 

In order to compare the results of this investigation 
with those published elsewhere on similar systems, convergence 
pressure was the only parameter to account for the effect 
of composition. Comparisons were made with the values of 
Kohn and Kurata (78) at Py, = 1800 psia and of Robinson et ay (139) 
and wesma (59) at Pyy = 1760 psia. Convergence pressure for 


this system appears to be the same as the critical pressure. 
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The methane K-values obtained in this investigation 
are slightly lower than the values reported by Kohn and 
Kereta at 140°F for the system methane - hydrogen sulphide. 
The hydrogen sulphide K-values cross at 1,300 psia. The 
values of Kohn and Kurata are higher for pressure less than 
1,300 psia but are lower than the values of this investigation 
for pressures above 1,300 psia. 

The methane K-values reported by Robinson et a ee 
compare well with the values obtained in this investigation. 
The curves of figure 11 cross at 1,300 psia. The values of 
Robinson et al are lower than the methane K-values presented 
here for pressure less than 1,300 psi and are higher above 
this pressure. 

The carbon dioxide K-values reported by Robinson et 
al are higher than the values obtained in this work but the 
difference between the two is not more than 0.3. The K-values 
equal each other for pressures less than 1,000 psia. The 
hydrogen sulphide K-values are always lower than the values 
obtained during this work. 

The methane K-values interpolated from the NGSMA 
Engineering Data Book for a convergence pressure of 1,760 


psia and at 120°F are much lower than the values of this 


investigation. 
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The carbon dioxide values compare well. The curves 
of figure 12 cross at 880 psia. The NGSMA values are 
slightly higher than the values obtained in this investi- 
gation for pressures above 880 psia and lower for pressures 
below this. 

The NGSMA values for hydrogen sulphide are always 


higher than the values presented here. 
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F. CONCLUSIONS 


(1) The two phase envelope of a quaternary natural 
gas system was established by measuring percent liquid present 
for different pressures at nine temperatures: 35°, 60°, 90°, 
E20R, LAO FuckSOShy LEO nel 702% candmieome. 

(2) Critical opalescence was observed and the 
phenomenon was compared with the observations of other workers. 

(3) The equilibrium ratios for hydrogen sulphide, 
carbon dioxide and methane were calculated at 120°F and compared 
with the values available from the literature on similar systems. 

(4) Due to the scanty data available on identical 
systems, the results of this investigation could only be 
compared with the data reported on H»9S-CHy and HzS-COs-CHy 
systems. It is not possible to make any general statements, 
but, comparison of the results serves to indicate that the 
convergence pressure is a potential and convenient parameter 
to account for the effects of composition on equilibrium 
ratios. 

(5) The existing equipment was modified and used for 
the first time to study the phase behaviour of a multicomponent 


system. It is suggested that the piston design be modified. 
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NOMENCLATURE 


fugacity of a component 

fugacity of a component in a mixture 

fugacity of the yth component in the pure liquid state 
fugacity of the jth component in the pure vapour state 
number of components in the system 

equilibrium ratio 

total pressure of the system 

vapour pressure of a component at the temperature of 
the equilibrium 

partial pressure of a component at the temperature and 
pressure of the equilibrium 

convergence pressure 

molal volume of vapour 

mole fraction in the liquid phase 


mole fraction in the vapour phase 


activity coefficient 


Subscripts 


a 


L 


Vi 


component i 
liquid phase 


vapour phase 
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Appendix I 


TABULATIONS OF EXPERIMENTAL RESULTS 
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Table 6 


TWO PHASE ENVELOPE FOR THE QUATERNARY NATURAL GAS SYSTEM 


Temperature Pressure Volume of the Volume of Percent 
ok psia System, ml Liquid ml Liguid 
35 246 0.00 0.00 

528 43.64 5.83 13.47 

778 ZS Oo 6.97 28.42 

1048 t2nos Vays 60.30 

L8u5 9257 8.47 88.50 

60 360 0.00 C.00 
590 49.47 5.43 UO. Sy 

TIS 20.54 

ES S055 6.53 21.00 

1000 40.63 

1024 7... 50 We 92 45.25 

1197 65.00 

1300 LB..82 9 Suk 78.76 

1365 10.87 9.16 84.26 

1540 9.42 8.99 95.45 

90 550 0.00 0.00 
750 48.65 4.78 9.82 

1011 25.44 6.90 2 iw 12 

1302 15.00 8.86 59.06 

1386 U3). US. 9.56 72.18 

1472 La fed LO, 21 80.07 

1552 12202 iO!. OM 90.76 

120 750 0.0 0.00 
S55 58.50 2.00 3.42 

950 45.88 3.44 H.50 

1200 26.91 6.18 22.96 

1477 16.94 8.86 52530 

155 14.88 9.79 65.80 
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Table 6 (Continued 


Temperature Pressure Volume of the Volume of Percent 
°F psia System, ml Liquid ml Liquid 
140 910 0.00 0.00 

1011 53.41 0.93 o/s 
372 256i: 4.03 IDES) 
1578 EO L55 5.98 30.60 
1748 15.46 siney/ 54.14 
1761 Critical Point 
150 1010 0.00 0.00 
1430 27250 2.20 8.06 
1789 LST) 1.65 9.29 
1640 19.94 Zuo 14.50 
1718 Poss Sa2s 17.14 
1335 0.00 0.00 
160 1135 48.97 0.05 OF. 16 
1318 36.70 Oe23 O63 
1821 187.93 0.35 E98 
1682 2dl0S 0.60 285 
1579 24.26 0.75 3.09 
1945 0.00 0.00 
170 1300 0.00 0.00 
1421 S57 7 6.15 0.42 
1615 PASI ROIAL Oe ay, 0.60 
aS)32 31.24 0320 0.64 
i7s2 24.76 O22 GO. Sl 
1930 0.00 0.00 
180 1560 SAO) 0.00 
1660 Zaye 0.10 0.40 
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Table 8 


VAPOR-LIQUID EQUILIBRIUM DATA FOR THE QUATERNARY 
NATURAL GAS SYSTEM AT 120°F 


Pressure Composition, Mole Fraction 
Sia Component Vapour Liquid 
926 H5S 0. 5562 0.9014 

CO2 0.1000 0.0498 
CH, 0.3300 0.0474 
No 0.0138 0.0014 
1030 H)S 0.5400 0.8481 
CO5 Or LLSO 0.0695 
CH, 0.3310 0.0808 
No 0.0140 0.0016 
1265 H5S 0.5429 0.8342 
CO2 0.1048 OL O7L5 
CHy Q2 3382 0.0924 
No 0.0141 0.0019 
1483 Hos 0.6753 0.7931 
CO5 0.0904 0.0782 
CHy 0.2288 0.1278 
N5 0.0055 0.0009 
1568 H5S ©: 7005 O21 7627 
CO» 0.0856 0.0844 
CHy 0.2023 6.1509 
No O 20006 0.0020 
1639 HS 0.6666 0.7543 
CO5 0.0929 0.0862 
CHy 0.2308 0.1582 
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Table 9 


EQUILIBRIUM RATIOS FOR THE QUATERNARY NATURAL 


Pressure 
psia 


926 


1030 


GAS SYSTEM AT 120°F 


Equilibrium Ratios 
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